Optimal Control of the Strong-Field Ionization of Silver Clusters in Helium Droplets 
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Optimal control techniques combined with femtosecond laser pulse shaping are applied to steer 
and enhance the strong-field induced emission of highly charged atomic ions from silver clusters 
embedded in helium nanodroplets. With light fields shaped in amplitude and phase we observe a 
substantial increase of the Ag 9+ yield for q > 10 when compared to bandwidth-limited and optimally 
stretched pulses. A remarkably simple double-pulse structure, containing a low-intensity prepulse 
and a stronger main pulse, turns out to produce the highest atomic charge states up to Ag 20+ . A 
negative chirp during the main pulse hints at dynamic frequency locking to the cluster plasmon. 
A numerical optimal control study on pure silver clusters with a nanoplasma model converges to a 
similar pulse structure and corroborates, that the optimal light field adapts to the resonant excitation 
of cluster surface plasmons for efficient ionization. 

PACS numbers: 36.40.Gk, 52.50.Jm 



I. INTRODUCTION 



Atomic clusters in intense laser pulses have received 
considerable attention since the mid-90's because of the 
potential for several promising applications, such as the 
controlled generation of highly charged atomic ions |Q, 
energetic electrons (2), x-rays ||, or even nuclear parti- 
cles Q. Besides technological interest, the absence of 
hidden dissipation channels makes clusters to a valu- 
able testing ground for exploring many-particle effects 
in strong-field laser-matter interactions || . 

A remarkable feature of clusters exposed to intense 
near-infrared (NIR) femtosecond laser pulses is the ex- 
tremely efficient absorption of radiation energy due to 
the coupling of the light field to a dense transient 
nanoplasma The resulting highly nonlinear cluster 
response has been studied theoretically by several groups 
and it is widely accepted that resonant collective electron 
heating by excitation of Mie surface plasmons is a key 
mechanism for high absorption as well as strong cluster 
ionization [M §1 Bl [r^, O. However, when considering 
NIR laser pulses, the energy of the plasmon resonance 
in metallic as well as in preionized rare-gas clusters is 
typically well above the laser photon energy of 1.5 eV in 
early stages of the interaction. For example, a plasmon 
energy of about 4.0 eV is found for silver clusters [[l2| in 
the ground state and may be further increased by laser- 
induced inner ionization. Noting the density dependence 
of the Mie-plasmon energy (w m i e ~ \fpi) 1 where p\ is the 
charge density of the ionic background) , resonant collec- 
tive heating requires a certain cluster expansion induced 
by ionization and non-resonant heating. Obviously, such 
sequence of ionization, expansion, and transient reso- 
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nance heating may be realized and resolved with very dif- 
ferent pulse shapes, e.g., by excitation with a long pulse 
or with multiple well-separated pulses, respectively [ p^| . 
Likewise, the interaction process can be sensitively tuned 
by modifying the laser parameters, e.g., by changing the 
pulse duration, envelope, or chirp. 

Indeed, several experiments have demonstrated that 
such simplified control techniques allow the enhancement 
of certain decay channels by optimizing the light field. 
In particular, cluster excitation with (i) stretched pulses 
having an optimal duration or with (ii) dual pulses having 
an optimal delay substantially support the gen eration of 
highly charged, energetic ions jl4| [ll| [r], [ItJ and fast 
electrons |l8], [l9|. At optimal pulse length, ion recoil 
energies can be further increased by (iii) introducing a 
weak additional shoulder in the leading edge of the laser 
pulse |2^] or by (iv) using a negative pulse chirp 

As a next logical step, optimal control techniques based 
on iterative feedback with genetic algorithms can be used 
to further explore and optimize the cluster excitation 
process. Considering the strong pulse-structure depen- 
dence obtained in previous studies (i-iv), laser-cluster in- 
teractions provide ideal grounds to test and apply con- 
trol schemes for steering the evolution of highly-excited 
many-particle system by a strong laser field. The first 
optimal control experiment on clusters in strong fields 
has been conducted by Zamith et al. (^2|, who used a 
liquid-crystal-based pulse shaper to maximize the abun- 
dance of highly charged Xe ions. The optimized light 
field contains a two-pulse sequence and produces ions up 
to Xe 23+ . The spectral phase of the optimized laser field, 
however, was not analyzed. 

In the present contribution we report on an optimal 
control experiment including a complete optimization 
and characterization of amplitude and phase of the laser 
field. As a prototype system wc study silver clusters em- 
bedded in helium nanodroplets, which have shown a high 




FIG. 1: Schematic view of the experiment: The top part visualizes the generation of the helium nanodroplets, the formation 
of the silver clusters inside the droplets by sequential atomic pickup in the oven, the laser excitation in the interaction region, 
and the time-of- flight analysis in the reflectron time-of- flight mass spectrometer (from left to right). The bottom part sketches 
the stages of pulse delivery from seed pulse generation in the oscillator, over pulse modification in the pulse shaper and final 
amplification, to pulse focusing into the molecular beam apparatus. In the optimal control experiment the yield of highly 
charged ions is maximized by applying an evolutionary algorithm scheme within a feedback loop (see arrows). 



pulse-structure sensitivity in previous experiments p5| . 
For this metallic system collective electron motion is ex- 
pected to dominate the optical absorption. In order to 
maximize the abundance of highly charged Ag 9+ ions a 
multi-parameter genetic algorithm procedure is used for 
tailoring the pulses. A sketch of the experimental proce- 
dure is displayed in Fig. [l| 

To evaluate the effect of the optimization, the re- 
sulting ion yields are compared to the spectra pro- 
duced with bandwidth-limited and dispersively stretched 
pulses. The optimal pulse shapes are characterized by 
frequency-resolved optical gating (FROG) Q and com- 
pared to results from a computational control experi- 
ment. The numerical optimization, based on a modified 
nanoplasma model varies the pulse envelope func- 
tion and is performed for a simplified system, i.e., pure 
silver clusters. 

The remainder of the text is organized as follows. Sec- 
tion H describes the experimental methods, i.e., the clus- 
ter generation, the laser system, and the details of pulse 
shaping and characterization. The results of the experi- 
mental and numerical closed-loop optimizations are pre- 
sented in Sec. The discussion and comparison of the 
data are subject of Sec. IV. Finally, conclusions are sum- 
marized in Sec. |v|. 



II. METHODS 



A. Experimental setup 



Silver clusters are produced by the helium droplet 
pickup technique p5[ , see Fig. |l|. Superfluid helium nan- 
odroplets are generated by means of a supersonic ex- 
pansion of highly purified helium gas (He 6.0) at a stag- 
nation pressure of 20 bars and a temperature of 8.9 K 
through a 5-fim nozzle yielding a mean droplet size of 10 6 
atoms. After differential pumping the molecular beam 
enters the pickup region where silver atoms are loaded 
into the droplet, forming clusters consisting of up to 150 
atoms |f26f . By changing the nozzle temperature and the 
metal vapor pressure in the pickup cell, the size of the 
droplet as well as the size of the metal core can be var- 
ied nearly independently over a broad range. We note 
that this type of cluster source shows an outstanding 
long-term stability, which is a central requirement for 
performing the optimization experiments. 

A ThSapphire laser system, consisting of the oscilla- 
tor, the pulse shaper, and the pulse amplifier, see Fig.[l], 
delivers pulses with energies of up to 2.5 mJ and pulse du- 
rations down to 35 fs (full width half maximum. FWHM) 
at 1 kHz repetition rate. Laser pulses enter the interac- 
tion region perpendicularly to the cluster beam and are 
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FIG. 2: Time-of-flight mass spectrum resulting from exci- 
tation of silver clusters in helium droplets by stretched 500 fs 
laser pulses at intensity 7 x 10 14 W/cm 2 (top panel: overview; 
bottom panel: enlarged section highlighting the range from 
Ag J+ to Ag 19+ ). The double peak structure of individual 
ion peaks results from the isotopic pattern of the silver atom 
( 107 Ag, 109 Ag). Due to the high recoil energies from the clus- 
ter Coulomb explosion the Ag ion features are broadened but 
can still be clearly resolved with the chosen detection tech- 
nique. In addition to Ag 9+ , also He + , He ++ , helium cluster 
ions, and Ag 9+ Hejv are observed. 



focused by a 30 cm lens (f# = 30) to intensities of up to 
10 16 W/cm 2 . Ions resulting from the laser-cluster inter- 
action are accelerated by a static electric field of 2 kV/ cm 
and analyzed by reflectron time-of-flight mass spectrom- 
etry p7[] . Fig.|| shows a typical mass spectrum resulting 
from excitation of the embedded metal complexes with 
500 fs pulses at intensity I = 7x 10 14 W/cm 2 . The pro- 
gression of multiply charged Ag 9+ ions with q up to 19 
can clearly be resolved with the applied detection tech- 
nique, see Fig.||(b). 



B. The Pulse Shaper 

The oscillator delivers seed pulses with a spectral band- 
width of 50 nm (FWHM) . Before amplification the pulses 
enter an acousto-optic programmable dispersive filter 
(AOPDF)(Dazzler, Fastlite) for pulse shaping, see Fig.[j]. 
The shaper modulates both the amplitude and the phase 



of the laser pulses simultaneously according to 

E out {u) oc S{oj)E in (u>), (1) 

with Ei n (uj) and E out (uj) the input and output laser fields 
in the frequency domain, and S(lo) the frequency re- 
sponse function of the filter. Arbitrary femtosecond laser 
pulse shapes can be generated with the appropriate S (uj) 
by controlling the acousto-optical signal. The frequency 
response function reads 

S(w) = A(u) exp[iB(w)], (2) 

where A(oj) and B(ui) are real functions describing the 
spectral amplitude and phase modulation, respectively. 
The AOPDF allows to realize a maximum group delay of 
3ps with 0.6 nm spectral resolution. For a more detailed 
description see Refs. |2^, p9[ . In our setup, A(lo) and 
B(uj) are defined by 60 parameters (30 in the range [0, 1] 
for A(oj) and 30 in the range [0, 27r] for B(cu)) equally 
spaced over the spectral range. 

In the control experiments these parameters are varied 
to optimize the pulse structure. We use an extended evo- 
lutionary algorithm scheme |30| to find the optimal laser 
pulse that maximizes the yield of highly charged atomic 
ions. For each pulse shape the mass spectrum is averaged 
over 6000 laser shots, whereas 50 individual laser pulse 
shapes enter each optimization cycle. The population 
of different pulse shapes is divided into separate sub- 
populations whose constituents are subject to crossover 
and random mutation after each generation. This pro- 
cedure has been proven to be more robust and efficient 
when compared to simple genetic schemes |3l| . The evo- 
lutionary process is iterated until the fitness value levels 
out. A full optimization experiment roughly takes 3-4 
hours. The final result of the optimization procedure, i.e., 
the laser pulse structure yielding the highest fitness value, 
is fully characterized in amplitude and phase by a home- 
build frequency-optical-gating (FROG) analyzer |2^] . 

III. RESULTS 

A. Experimental findings 

A comparison of the ion spectra resulting from excita- 
tion of the embedded silver clusters with (i) bandwidth- 
limited, (ii) optimally stretched, and (iii) fully shaped 
laser pulses is given in Fig.||. The data in Fig.||(a) shows 
a spectrum for cluster excitations with 50 fs pulses (i) at 
peak intensity Iq = 10 16 W/cm. In this case the ion sig- 
nal is dominated by Hcjv ions and weakly charged species 
of residual gas molecules (e.g., C 2+ ). Zooming in, how- 
ever, reveals the generation of highly charged silver ions 
up to Ag 17+ from the clusters with low yields, e.g., with 
about 1.5 x 10 -3 ions per laser shot for Ag 17+ . From 
previous studies it is well-known that longer pulses en- 
hance the ion yield, see e.g. (l4). We utilize the AOPDF 
to stretch the pulses (ii) by introducing a linear temporal 
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FIG. 3: Mass spectra with highly charged atomic ions ob- 
tained when exposing embedded silver clusters to intense laser 
pulses of different pulse shapes: (a) 50 fs pulses; (b) 500 fs lin- 
early down-chirped pulses (cf. Fig. ^) ; (c) laser pulses fully 
optimized in amplitude and phase within the control exper- 
iment. Note, that the scales are directly comparable. Ob- 
viously the yield of Ag 9+ is substantially enhanced with the 
tailored pulses. The inset in (c) shows the evolution of the 
average fitness value in the optimization experiment, i.e., the 
total yield of Ag 9+ (q=14-17) (open squares) compared to the 
reference measurements (50 fs: full square; 500 fs pulses: full 
triangle) . 

down-chirp on the frequency component. This results in 
a substantial increase in the yield of multiply-charged 
species from the clusters while ion contributions from 
background gas decline. At a pulse width of about 500 fs 
(I = 10 15 W/ cm, -2 ), see Fig.||(b), we observe the max- 
imal signal of Ag 9+ , indicating optimal ionization condi- 
tions with a single pulse. Ag 19+ shows up in the spectrum 
as the highest charge state. The two single-pulse mea- 
surements (i, ii) serve as references for the optimization 
experiments. 

For the full optimization procedure (iii) the yield of 
Ag 9 + (q = 14—17) is taken as fitness value for the genetic 
algorithm. The resulting final mass spectrum is given in 
Fig. H(c). The evolution of the fitness value (see inset) 
shows that the pulse structure quickly adapts to the de- 
tails of the cluster dynamics. The fitness value exceeds 
the value for the optimally stretched single pulse (solid 
triangle) already after a few iterations. Subsequently, the 
fitness value begins to level out and saturates at about 
twice the stretched pulse value. Hence, the high-g ion 
yield with the fully shaped light field is enhanced by 
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FIG. 4: Pulse shape (solid) and evolution of the carrier wave- 
length (dashed) of the laser fields corresponding to the three 
excitation scenarios in Fig. H (a) 50 fs pulse; (b) 500 fs linearly 
down-chirped pulse; (c) fully optimized light field. The tai- 
lored laser field leading to the enhanced generation of highly 
charged ions exhibits a double pulse structure, comprised of 
a weak, slowly rising prepulse and a stronger main pulse sep- 
arated by 140 fs. 



more than a factor of two when compared to the opti- 
mal single pulse, see Fig.||(b) and (c). The increase is 
even more pronounced (factor 17) when comparing to 
the spectra for the shortest and thus most intense pulse, 
cf. Fig. ||(a). In addition, we also observe an increase 
in the highest charge state reaching q max = 20 with the 
fully optimized pulse shape. 

As a first step towards unraveling the underlying dy- 
namics, the pulses from the three excitation scenarios in 
Fig. ^|(a)-(c) have been characterized in amplitude and 
phase, see Fig.|4](a)-(c). For convenience, the temporal 
intensity is given (solid) together with the phase evo- 
lution in terms of the instantaneous carrier wavelength 
(dashed) . which is derived from the instantaneous angu- 
lar frequency u)(t) = luq — </>(£), with u>o being the central 
angular frequency and <p(t) the temporal phase of the 
pulse g|. 

Whereas the flat evolution of the instantaneous car- 
rier wavelength in Fig.|^(a) underlines excitation with 
a nearly bandwidth-limited pulse, the linear wavelength 
shift in Fig. [|(b) reflects the dispersive pulse stretching. 
The analysis of the optimized laser field, see Fig. ||(c), 
shows a simple double-pulse intensity distribution (solid 
line). A prepulse containing about 15% of the pulse en- 
ergy excites the clusters prior to the main pulse. The 



5 



pulse intensity drops significantly between the subpulscs, 
indicating a two-step excitation process to induce the ef- 
fective generation of high charge states. This signature of 
the pulse is found to be rather robust with respect to the 
chosen pulse energy, i.e., the pulse delay as well as the 
intensity ratio may change while the general shape is pre- 
served. Note that the optimization procedure introduces 
a substantial red-shift of the instantaneous carrier wave- 
length from 790 nm to 810 nm during the main pulse, see 
dashed curve in Fig. |](c). 



B. Computational results 




To assist the experimental study, we apply a theoret- 
ical approach to find laser pulse shapes that maximize 
the yield of highly charged ions in the laser-cluster in- 
teraction. For this purpose, a genetic algorithm is used 
and adapted to the nanoplasma approach of Ditmire . 
The model contains the essential physical processes like 
ionization, heating and expansion. Briefly, the clusters 
are assumed to have uniform temperature and density 
profiles and are initialized as neutral spheres. To model 
the nanoplasma inner ionization dynamics, tunneling as 
well as electron impact ionization are taken into account. 
Heating via inverse bremsstrahlung, being the dominant 
process for laser energy capture, is described with effec- 
tive rates for collisional absorption. The cluster expan- 
sion and the resulting cooling of the system are calculated 
using hydrodynamic equations. In our version of the 
nanoplasma model, plasma heating is treated by a gen- 
eralized quantum statistical expression for electron-ion 
collisions including resonant collective absorption p2| . 
In addition, damping due to electron collisions with the 
cluster boundary, i.e., surface friction, is incorporated by 
including a corresponding collision term |?3|, Fur- 
thermore, the lowering of the ionization energies due to 
screening is taken into account Q . 

To calculate the optimal light field for producing a 
maximum yield of specific ion species a genetic feedback 
algorithm is applied. For the optimization procedure, a 
pulse intensity shape I(t) according to the Ansatz 



I(t) = ^2 A n exp 



n=0 



41n2(t-t r 



(3) 



is used, with A n: t n , and a n being the weights, time 
offsets, and temporal width parameters (FWHM) of the 
individual terms. Note that a fixed carrier wavelength 
of 810 nm is used in the calculations, i.e., the phase evo- 
lution is not optimized. The total pulse fluence is kept 
constant for each chosen pulse-shape and corresponds to 
a 810 nm Gaussian pulse of 35 fs (FWHM) with peak in- 
tensity 1 x 10 16 W/cm 2 . Further, a lower limit of 35 fs is 
used for the width parameters cr n . To calculate the best 
pulse shape, 20 different pools, each with a population of 
40 species, are used in the optimization procedure. 



FIG. 5: Calculated temporal intensity evolution of the opti- 
mal light field for producing a maximum yield of Ag 10+ . The 
intensity envelope is in qualitative accordance with the exper- 
imental result, cf. Fig. tA. The total pulse fluence corresponds 
to a 35 fs pulse with peak intensity 1 x 10 16 W/cm 2 . 



We perform the calculations as close as possible to 
the experimental conditions. Silver clusters with an ini- 
tial radius of 2.5 nm (4500 atoms) are considered since 
smaller particles cannot be treated adequately within the 
nanoplasma model. The calculated optimal pulse shape 
giving the maximum yield of Ag 10+ is shown in Fig. 0. 
This ion species was chosen because higher ionization lev- 
els were only weakly populated. Here, further improve- 
ments of the model would be necessary. For a discussion 
of additional processes to be included, like the enhance- 
ment of electron impact ionization by local cluster fields 
and the influence of the experimental setup, especially 
the frustration of recombination by the ion extraction 
electric fields after the laser pulse, see Rcf. p5| . 

As references to assess the results of the genetic algo- 
rithm, calculations were also performed for a single 500 fs 
pulse and for a dual-pulse excitation scheme. Whereas 
the single pulse mainly produces Ag 7+ and Ag 8+ , Ag 8+ 
and Ag 9+ are the dominant contributions in the ion spec- 
tra for dual pulse excitation. With the light field from the 
computational control experiment, see Fig. ||, the Ag 10+ 
yield is enhanced by a factor of four over the result for 
dual-pulse excitation and by more than a factor of eight 
over the result for the 500 fs single pulse. The simulations 
predict an optimized light field having a low-intensity 
prepulse followed by a more intense main pulse, in good 
agreement with the experimental control study. Even 
the intensity ratio between the two subpulses is nearly 
reproduced by the calculation. 



IV. DISCUSSION 

The applied pulse shaping technique combined with a 
genetic algorithm has been shown to allow the enhanced 
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production of highly charged ions from strong-field laser- 
cluster interactions. Compared to simpler excitation 
schemes with bandwidth-limited or stretched pulses, the 
ion yields could be substantially increased. The laser field 
tailored in amplitude and phase thus offers an alternative 
and more efficient ionization pathway. 

We find that a sequence of a weak prcpulsc followed by 
a much stronger main pulse after about 140 fs leads to the 
highest ionization under the experimental boundary con- 
ditions. Interestingly, the bandwidth-limited pulse from 
the reference scenario and the main feature of the op- 
timized light field are similar in intensity and duration. 
This underlines the decisive effect of the low-intensity 
prepulse (f0 15 W/cm 2 ) for enhancing the ionization effi- 
ciency by up to an order of magnitude. The optimized 
light field also outperforms optimally stretched single 
pulses (500 fs) by a factor of two in the ion yield. 

The envelope of the tailored light field and comparison 
with the results for a bandwidth-limited pulse suggest 
that efficient energy transfer and strong ionization take 
place within the stronger main pulse. These findings can 
be related to previous dual pulse studies on Ag^r in He 
droplets , where a specific optical delay was found to 
maximize the highly charged ion yield for excitation with 
a pair of identical pulses. The optimization presented 
here extends these experiments. The control study not 
only underlines the importance of a certain pulse delay 
that matches the cluster expansion required for most effi- 
cient resonant collective heating, but also highlights the 
impact of the pulse intensity ratio. The genetic opti- 
mization shows that a less intense prepulse is favorable 
for preionization of the particle in order to launch the 
cluster expansion. Near resonance, the pulse adapts to 
maximum heating and strongest charging. This general 
scenario is corroborated by the fact, that the observed 
characteristic pulse structure is quite robust with respect 
to changes of the pulse energy 

The observed pulse shape shows similarities with the 
result of an earlier control experiment by Zamith et 
al. HH on Xc N (N w 1.6 x 10 4 ). The authors also found 
an asymmetric dual pulse with a time separation of f 46 fs. 
Although the pulse envelope indicates a similar scenario, 
the actual target conditions are quite different. In silver 
clusters the delocalized valence electrons already form a 
metallic state, whereas a rare-gas system requires preion- 
ization for the nanoplasma formation. 

Comparing the outcome of the present experimental 
optimization with the numerical one shows similar pulse 
structures. A detailed analysis of the simulation run for 
the numerically optimized light field corroborates that 
the main energy transfer into the system takes place upon 
the passage of the plasmon resonance in the trailing main 
pulse. Here, the major fraction of the energy capture oc- 
curs leading to strong heating and finally to the formation 
of a large number of highly charged ions. The compu- 
tational results thus support that the optimal light field 
" tries" to utilize most efficient plasmon enhancement ion- 
ization Q. 

Interestingly, the analysis of the instantaneous phase 



reveals a red-shift of the central carrier wavelength 
(down-chirp) during the main laser pulse, cf. Fig. ^(c). 
This behavior can be explained with dynamic frequency 
locking of the laser field to the plasmon. As the clus- 
ter expands, a down-chirped laser field allows extended 
frequency matching with the collective mode. For this 
to occur, however, the increase in the ion density due to 
cluster inner ionization has to be overcompensated by a 
sufficiently rapid cluster expansion. 

Finally, we like to discuss the impact of the helium 
nanomatrix. Whereas helium is almost transparent at 
infrared wavelengths for low intensities, substantial ion- 
ization and heating may occur under the strong laser 
field. In our calculations the helium is not taken into 
account up to now due to the restricted dimensional ap- 
proach of the nanoplasma model. Very likely, the sur- 
rounding helium droplet influences both the timing as 
well as the final ionization state of the cluster ions. Fur- 
ther improvements of the model for resolving the radial 
structure of the system are underway. As has been re- 
ported recently by Mikaberidze et al. 11 for embedded 
xenon clusters (XeiooHeiooo), the presence of the helium 
nanomatrix results in a large contribution to the total 
energy absorption in a strong laser field. The impact of 
the helium matrix on the final charge state distribution 
of the cluster constituents as probed in our experiment, 
however, remains a challenging problem for future inves- 
tigations. 



V. CONCLUSION 

Amplitude- and phase-modulated pulses have been 
used to maximize the highly charged ion yield from silver 
clusters embedded in helium nanodroplets. A substantial 
increase is found when comparing to bandwidth-limited 
and optimally stretched pulses. The laser field resulting 
from the control experiment shows a two-pulse structure 
which is in fair agreement with the result of a computa- 
tional optimization experiment based on the nanoplasma 
model. The shape and phase evolution of the optimized 
laser field indicate that the tailoring adapts to take ad- 
vantage of the plasmon enhanced response. This proof 
of principle experiment on the metallic prototype system 
thus reveals a unique signature from collective electron 
excitations and their importance for controlling strong- 
field laser-matter interactions. 



Acknowledgments 

Main parts of the helium droplet machine have 
been provided by J. P. Toennies and his group at the 
MPI Gottingen. Financial support by the Deutsche 
Forschungsgcmcinschaft within the Sondcrforschungs- 
bereich SFB 652 is gratefully acknowledged. 



7 



[1] E. M. Snyder, S. A. Buzza, and A. W. Castleman Jr., 

Phys. Rev Lett. 77, 3347 (1996). 
[2] E. Springate, S. A. Aseyev, S. Zamith, and M. J. J. 

Vrakking, Phys. Rev. A 68, 053201 (2003). 
[3] A. McPherson, B.D. Thompson, A.B. Borisov, K. Boyer, 

and C.K. Rhodes, Nature 370, 631 (1994). 
[4] T. Ditmire, J. Zweiback, V.P. Yanovsky, T.E. Cowan, 

G. Hays, and K.B. Wharton, Nature 398, 489 (1999). 
[5] U. Saalmann, C. Siedschlag, and J. M. Rost, J. Phys. B 

39, R39 (2006). 
[6] T. Ditmire, R. A. Smith, J. W. G. Tisch, and M. H. R. 

Hutchinson, Phys. Rev. Lett. 78, 3121 (1997). 
[7] T. Ditmire, T. Donnelly, A. M. Rubenchik, R. W. Fal- 
cone, and M. D. Perry, Phys. Rev. A 53, 3379 (1996). 
[8] H.M. Milchberg, S.J. McNaught, and E. Parra, Phys. 

Rev. E 64, 056402 (2001). 
[9] M. Kundu and D. Bauer, Phys. Plas. 15, 033303 (2008). 
[10] U. Saalmann and J.-M. Rost, Phys. Rev. Lett. 91, 223401 

(2003). 

[11] A. Mikaberidze, U. Saalmann, and J. M. Rost, Phys. 

Rev. A 77, 041201(R) (2008). 
[12] J. Tiggesbaumker, L. Koller, K.-H. Meiwes-Broer, and 

A. Liebsch, Phys. Rev. A 48, R1749 (1993). 
[13] J. Zweiback, T. Ditmire, and M. D. Perry, Phys. Rev. A 

59, R3166 (1999). 
[14] L. Koller, M. Schumacher, J. Kohn, S. Teuber, 

J. Tiggesbaumker, and K.-H. Meiwes-Broer. Phys. Rev. 

Lett. 82, 3783 (1999). 
[15] T. Doppner, Th. Fennel, Th. Diederich, 

J. Tiggesbaumker, and K.-H. Meiwes-Broer, Phys. 

Rev. Lett. 94, 013401 (2005). 
[16] V. Kumarappan, M. Krishnamurthy, and D. Mathur, 

Phys. Rev. A 66, 033203 (2002). 
[17] A. Heidenreich, I. Last, and J. Jortner, Las. Phys. 17, 

608 (2007). 

[18] Th. Fennel, T. Doppner, J. Passig, Ch. Schaal, 
J. Tiggesbaumker, and K.-H. Meiwes-Broer, Phys. Rev. 
Lett. 98, 143401 (2007). 

[19] V. Kumarappan, M. Krishnamurthy, and D. Mathur, 



Phys. Rev. A 67, 043204 (2003). 
[20] K. Mendham, J. Tisch, M. Mason, N. Hay, and J. Maran- 

gos, Opt. Express 11, 1357 (2003). 
[21] Yuji Fukuda, Koichi Yamakawa, Yutaka Akahane, 

Makoto Aoyama, Norihiro Inoue, Hideki Ueda, and Ya- 

suaki Kishimoto. Phys. Rev. A 67, 061201(R) (2003). 
[22] S. Zamith, T. Martchenko, Y. Ni, S.A. Aseyev, H.G. 

Muller, and M.J.J. Vrakking, Phys. Rev. A 70, 

011201(R) (2004). 
[23] R. Trebino, Frequency -Resolved Optical Gating: The 

Measurement of Ultrashort Laser Pulses. Kluwer Aca- 
demic, London (2000). 
[24] P. Hilse, M. Moll, M. Schlanges, and Th. Bornath, Las. 

Phys. 19, 428 (2009). 
[25] A. Bartelt, J.D. Close, F. Federmann, N. Quaas, and J.- 

P. Toennies, Phys. Rev. Lett. 77, 3525 (1996). 
[26] Th. Diederich, J. Tiggesbaumker, and K.-H. Meiwes- 
Broer, J. Chem. Phys. 116, 3263 (2002). 
[27] T. Doppner, Th. Diederich, A. Przystawik, N.X. Truong, 

Th. Fennel, J. Tiggesbaumker, and K.-H. Meiwes-Broer. 

Phys. Chem. Chem. Phys. 9, 4639 (2007). 
[28] F. Verluise, V. Laude, Z. Cheng, Ch. Spielmann, and 

P. Tournois. Opt. Lett. 25, 575 (2000). 
[29] DAZZLER, http://www.fastlite.com. 
[30] H. Pohlheim, Evolutionare Algorithmen: Verfahren, Op- 

eratoren und Hinweise fur die Praxis. Springer, Berlin 

(2000). 

[31] David E. Goldberg, Genetic algorithms in search, opti- 
mization, and machine learning. Addison- Wesley Profes- 
sional, New York (1989). 

[32] Th. Bornath, P. Hilse, and M. Schlanges, Contrib. 
Plasma Phys. 17, 591 (2007). 

[33] F. Megi, M. Belkacem, M. A. Bouchene, E. Suraud, and 
G. Zwicknagel, J. Phys. B 36, 273 (2003). 

[34] S. Micheau, H. Jouin, and B. Pons, Phys. Rev. A 77, 
053201 (2008). 

[35] Th. Fennel, L. Ramunno, and Th. Brabec, Phys. Rev. 
Lett. 99, 233401 (2007). 



